THE TILT OF THE FUNDAMENTAL PLANE: 
THREE-QUARTERS STRUCTURAL NONHOMOLOGY, 
ONE-QUARTER STELLAR POPULATION 
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ABSTRACT 

The variation of the mass-to-hght ratios M/L of early type galaxies as function of their lumi- 
nosities L is investigated. It is shown that the tilt j3 = 0.27 (in the B-band) of the fundamental 
plane relation M/L (x can be understood as a combination of two effects: about one-quarter 
(i.e. Af3 = 0.07) is a result of systematic variations of the stellar population properties with in- 
creasing luminosity. The remaining three-quarters (i.e. A(3 = 0.2) can be completely attributed 
to nonhomology effects that lead to a systematic change of the surface brightness profiles with 
increasing luminosity. Consequently, the observed tilt in the K-band (/3 = 0.17) where stellar 
population effects are negligible, is explained by nonhomology effects alone. After correcting for 
nonhomology, the mean value of the mass-to-light ratio of elliptical galaxies (M/L^) is 7.1±2.8 
(Icr scatter). 



Subject headings: galaxies: clusters: general — galaxies: interactions 
photometry — galaxies: fundamental parameters 



.laxies: structure — galaxies: 



1. Introduction 

Traditionally, the global properties of elliptical 
galaxies have been characterized by three observ- 
ables: the central velocity dispersion co, the effec- 
tive or half-light radius re, and the mean surface 
brightness within the effective radius /g. It is gen- 
erally assumed that elliptical galaxies are relaxed 
particle systems. From a theoretical point of view, 
the virial theorem then predicts a relationship be- 
tween the gravitational mass, the velocity disper- 
sion and the radius. This can be written as a rela- 
tionship between /e , ctq and r^ . These observables 
should then be constrained to a two-dimensional 
plane in the above three-space. Elliptical galax- 
ies are indeed found to lie along a plane; how- 
ever, the observed plane (popularly known as the 
fundamental plane [FP]; Djorgovski & Davis 1987, 
Faber et al. 1987, Dressier et al. 1987) deviates 
from the prediction using a constant M/L. This 



deviation is known as the "tilt" of the FP. Since 
the above prediction is based on the assumption 
that the mass-to-light ratio and the global struc- 
ture of the elliptical galaxies are not dependent on 
the luminosity, the explanations for the tilt have 
been attributed to deviations from one or both of 
these hypotheses. 

Following on that, Faber et al. (1987) sug- 
gested that the M / L values increase systemati- 
cally with increasing luminosity. Using this inter- 
pretation, the observed averaged value of the FP 
tilt over a large group of galaxy samples in the B- 
band is written as M/L cx i'^, with /3=0.27±0.08 
(Prugniel & Simien 1996). Since brighter ellip- 
ticals are systematically redder, at least part of 
the tilt must be understood as stellar population 
variations along the different luminosities of ellip- 
tical galaxies. However, the presence of a "re- 
maining" tilt (/3=0.17±0.01) in the K-band FP 
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(Pahrc, Djorgovski & de Carvalho 1998), where 
population variations are weak, implies that this 
solution alone cannot explain completely the ob- 
served slope. In fact, the variation in the prop- 
erties of the stellar population as deduced from 
broad band colors and the Mg2 index explain only 
one-third of the slope, i.e., A/3 0.10 in the B- 
band (Tinsley 1978, Dressier et al. 1987, Prugniel 
& Simien 1996). We show this more explicitly in 
Figure la, where the luminosity dependence of 
estimated stellar M/L ratios in B (open circles) 
and JsT-bands (filled circles) is plotted for a sample 
of 911 morphologically selected early- type galaxies 
from the Sloan Digital Sky Survey (SDSS) Early 
Data Release and the Two Mass All Sky Surveys 
(2MASS) (BeU et al. 2003). The dependence of 
stellar M/L ratios on galaxy luminosity is weak 
in the B-band, in line with previous studies, and 
is negligible in the /C-band. Taken together, the 
above is consistent with the conclusion that the 
dominant part of the tilt of the FP (A/? ~ 0.2) is 
not a population effect. 

Changes in the mixture of dark-to- visible mat- 
ter as function of luminosity would be a plausible 
explanation for the remaining tilt. However, ob- 
servations (Romanowsky et al. 2003) demonstrate 
that the amount of dark matter within an effec- 
tive radius is in general negligibly small. Further- 
more, systematic variations in kinematical struc- 
ture (van der Marel 1991, Bender, Burstein and 
Faber 1992, hereafter BBF) and radial anisotropy 
(Ciotti, Lanzoni & Renzini 1996, Ciotti & Lanzoni 
1997) cannot produce the observed tilt. The role 
played by rotation is also small (Busarello et al. 
1997). 

In this Letter, we argue that the remaining tilt 
of the FP is due to a systematic variation of the 
structural (and its associated dynamical) nonho- 
mology of the elliptical galaxies as function of their 
luminosity. Following an analogy with the FP 
of the globular clusters, Djorgovski (1995) sug- 
gested a systematic change in the galaxy concen- 
tration as a viable hypothesis. Hjorth & Mad- 
sen (1995) explained part of the tilt of the FP 
(A/3 = 0.12) using a galaxy model with a variable 
central gravitational potential correlated with the 
luminosity. Following a similar idea but using a 
model directly related with the observations (the 
Sersic (1968) r^/" law), Ciotti et al. (1996) sug- 
gested that a suitable variation (by a factor of ^ 



2 -4) of the shape parameter n (i.e. a different de- 
gree of concentration) along the FP would be able 
to explain the tilt. However, Prugniel & Simien 
(1997) (and see also Graham & CoUess (1997) in 
the V-band) showed that the structural nonho- 
mology contributes only A/3 ~ 0.06 — 0.08 and to 
fully account for the tilt of the FP they needed 
to include some contributions from stellar popula- 
tion and rotational support. These analysis could 
be affected by a limited sample which had a small 
range in Sersic indices 2<n<5 which is far from 
providing the full range of nonhomology observed 
by other authors with l^n<10 in the same abso- 
lute magnitude range -22<Mb<-15 (see e.g., Gra- 
ham & Guzman 2003). 

Taking into account the full range of structural 
nonhomology observed, we quantify the contribu- 
tion to the tilt from the structural nonhomology 
in the B band. We show that the structural non- 
homology is able to explain «3/4 of the tilt in the 
B-band (i.e. A/3=0.2), which, in combination with 
the tilt expected from population variations (i.e. 
A/3~0.1) can explain completely the observed tilt 
of the FP. 

2. The scalcir virial theorem and the FP 
tilt 

We use the identity 

L = cihrl (1) 

and the scalar virial theorem for a stationary stel- 
lar system 

GM = c-ialr, (2) 

with L being the luminosity, G the gravitational 
constant, M the mass and c\ and C2 structure 

terms. Any stellar system in virial equilibrium is 
then expected to follow the identity: 

r, = G-\c^c{^){MlL)-^oll-\ (3) 

With the definition given in the Introduction for 
/e, ci = 27r. The FP (Dressier et al. 1987; Djor- 
govski & Davis 1987) in the B-band can be fitted 

r-e « [airi-^^'. (4) 

This relation is a consequence of the variation 
of c^'^iM/L) with L} If, following Faber et al. 

In the more general case, Cj -^(M/L)=Ll'^l2, 7 has been 
found to be compatible with (7 = —0.01 ± 0.13; Prugniel 
& Simien 1996) and in what follows we will assume 7 = 0. 
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(1987); C2 is constant elliptical galaxies would 
be structurally homologous systems and the tilt 
would have to be explained as a systematic varia- 
tion of mass-to -light (M/L) ratio with luminosity: 
M/L oc Lf^, with /3=0.27±0.08. 

As discussed in the Introduction it is unlikely 
that this tilt can be explained by stellar popula- 
tion effects alone. One needs to combine two well 
known observations: the systematic variation of 
stellar populations with luminosity and the change 
of observed structural properties along the early- 
type sequence from dwarf to giant ellipticals. Ide- 
ally both effects in combination would be able to 
account for the full tilt of the FP. 

3. The relation between C2 and n 

It has been shown (see e.g. Trujillo, Graham & 

Caon 2001 and references therein) that elliptical 
galaxies do not form a homologous structural fam- 
ily and that the luminosity-dependent departures 
from the r^/^ law can be described by the r-'^/" 
Sersic model. The nonhomology is also reflected in 
the strong correlations between the shape param- 
eter n and photometric-independent galaxy prop- 
erties as, for example, the central velocity disper- 
sion (Graham, Trujillo & Caon, 2001). 

In Fig. lb, we show the relation between the 
shape index n and the absolute B-band (model- 
independent) magnitude for 200 elliptical galax- 
ies. The galaxies used in this plot correspond 
to ellipticals from the Virgo, Fornax and Coma 
Clusters (Caon et al. 1990; Caon, Capaccioli & 
D'Onofrio 1994; Binggch & Jerjen 1998; Gutierrez 
et al. 2004). Those galaxies classified as SO were 
removed from our sample to avoid possible miss- 
measurements of the index n due to the disk com- 
ponent of these galaxies. Error estimates for n are 
found to have a typical uncertainty of 25% (Caon, 
Capaccioli & D'Onofrio 1993). The high statisti- 
cal significance of this correlation has been studied 
in previous papers (see, e.g. Graham et al. 2001). 
To check that the robustness of the above relation 
is not affected by the large number of faint galax- 
ies, we have evaluated the Spearman correlation 
coefficient for those galaxies brighter than Mb=- 
16.5 (Ho=70 km s"! Mpc"!) and we find rs=0.72. 

Using the Sersic model, it is possible to relate 
the structural parameter C2 with n, leading also to 
a relation between n and L, or between C2 and L. 



To estimate C2 for a Sersic model we need to eval- 
uate the velocity dispersion profile ^^(r). In the 
most simple case this can be done by assuming a 
spherical, nonrotating, isotropic i^/"^ model and 
then evaluating Gr{r) using the Jeans equation 
(see, e.g., Binney & Tremaine 1987). However, 
the observed central velocity dispersion (To that 
enters in Eq. 2 does not correspond to crr(O), but 
rather to the observed projected velocity disper- 
sion (Tp (i?) , luminosity averaged over the aperture 
used for the spectrographic observations crap{Rap)- 
Consequently, to evaluate properly C2, we need 
to switch from Uri^) to (Jp{R) and from this to 
(^ap{Rap) (see a detailed explanation e.g. in Sec. 
2.2 of Ciotti et al. 1996). 

To ilhistrate how C2 changes depending on 
which estimate of the velocity dispersion is used 
we show in Fig. 2 the relation between C2 and n 
derived from ar{re/8), (Tp(re/8) and aap{re/&)- 
We show C2 for r = re/8 because this is one 
of the most usual aperture radii to measure cto 
(J0rgensen, Pranx & Kj^rgaard 1993). Our value 
of C2 is in good agreement with the estimates 
of this quantity by other authors (e.g. Prugniel 
& Simien (1997), Bertin, Ciotti & Del Principe, 
2002). 

Since the sizes of galaxies differ from one an- 
other and because galaxies are observed at differ- 
ent distances, the estimation of the central velocity 
dispersion using a fixed angular aperture samples 
different fractions of the total light (or the effective 
radii). Consequently, when dealing with samples 
that extend over a large range of sizes, it is not 
the best approximation to estimate C2 from the 
observations by using a fixed aperture related to 
the effective radius (e.g. re/8). Instead one should 
use typical angular apertures (e.g. l."6 or 2. "2). 

In order to evaluate the influence of nonhomol- 
ogy in the FP, we have selected from the sample 
of elliptical galaxies presented above those galax- 
ies that have a measured central velocity disper- 
sion. The velocity dispersions are obtained from 
Hypercat. This leaves us with a total of 45 galax- 
ies ranging from -15 to -22 in the B band. In 
Fig. 2 we show the C2 values for the galaxies of 
this subsample estimated using two different fixed 
realistic apertures l."6 and 2. "2. As we can see 
in this figure, the effect of changing the aperture 
from l."6 to 2. "2 is not dramatic. In what follows 
we will use the mean value of these two measure- 
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mcnts C2(l."6) and C2(2."2) as an estimate of the 
true value of C2 for our observed galaxies. 

4. The tilt of the FP 

Fig. 3 shows the relation between M/L and the 
luminosity L for the selected galaxies, taking into 
account the variation in C2 as a result of variations 
in shape parameters n. For comparison, in the 
same figure we also show the result for a constant 
value for C2, i.e. if structural homology is assumed 
to hold. The value we have used for C2=4.86 in 
the assumption of homology is taken to match the 
value presented in BBF for a King model with 
I't/i'c ~100. This value is quite close to the value 
expected for a dc Vaucouleurs law (n=4; C2=4.6). 

Fig. 3 shows that the relationship between the 
M/ L and L is much flatter if the nonhomology is 
accounted for. To quantify this change we have 
performed a minimum fit to the data distri- 
butions using the relation {M/L) = aL^ . The 
values are (3 = 0.10(±0.04) if nonhomology is ad- 
dressed and P = 0.29(±0.04) for the homologous 
case. The exponent /3 for the homologous case 
in the B-band is in excellent agreement with the 
values reported previously. On the other hand, 
we can see that the nonhomology part of the tilt 
accounts for A/? ^ 0.2. We have also checked if 
the incompleteness is playing a role at the fainter 
luminosities. Our sample is 90% complete down 
Mb=-17. If we remove from our analysis all the 
galaxies below this limit we get the following ex- 
ponents: /3 = 0.27(±0.04) for the homologous case 
and (3 = 0.06(±0.04) in the nonhomologous case. 
Again, we find that the nonhomology part of the 
tilt accounts for A/3 ~ 0.2. 

5. Discussion 

It is encouraging that the tilt (3 ~ 0.1 of the 
FP after correction for to nonhomology is in very 
good agreement with the value expected due to 
systematic variations of stellar population proper- 
ties of elliptical galaxies (A/3 ~ 0.07). The ob- 
served tilt of /3 = 0.29(±0.04) in the B band FP 
therefore can be understood as a combination of 
three-quarters structural nonhomology and one- 
quarter stellar population effects. The tilt of the 
K-band FP where the change in the stellar popu- 
lation is negligible is approximately /3 = 0.2. This 
is indeed in good agreement with our expectation 



if the tilt were due only to the nonhomology and if 
color gradients inside galaxies do not change their 
Sersic indices. Indeed, radial V-K gradients in 
the surface brightness profiles of ellipticals are on 
average only 0.16 mag dex~^ (Peletier, Valentijn 
& Jameson 1990) which is small compared to the 
variations in surface brightness within an effective 
radius. The homologous FP of our sample has a 
scatter in log(M/L) of order 25% which is consis- 
tent with previous results (BBF; Pahre, Djorgov- 
ski & Carvalho 1998). The scatter is not signifi- 
cantly reduced in the non-homologous FP which 
indicates that it is mainly a result of population 
effects and not of structural variations within a 
given luminosity bin. 

Our results rule out that even very massive el- 
lipticals are strongly dominated by dark matter 
in their inner regions, which is in agreement with 
Rix et al. (1997) and Gerhard et al. (2001). Af- 
ter correcting for nonhomology, the mean value 
of the mass-to- light ratio (M/L^) of early- type 
galaxies in our luminosity range is 7.1±2.8 (Icr 
scatter). This value is in good agreement with 
the value expected for an old stellar populations 
of M/Ls=7.8±2.7 (Gerhard et al 2001). 

We thank Nicola Caon and Helmut Jerjen for 
providing us with an (updated) electronic copy of 
the Virgo Cluster Catalogs and Hans- Walter Rix, 
Marijn Franx, Alister W. Graham and Ralf Ben- 
der for interesting discussions. Some data used 
in this work were been obtained from HyperLeda 
database. 
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Fig. 1. — {a) Luminosity dependence in stellar 
M/L ratios in the -B-band (open circles) and the 
K-h&nd (filled circles) from 911 morphologically 
selected early-type galaxies (14.5 < r < 16.5; K < 
13.6; re > 2") from a combined SDSS/2MASS 
galaxy catalog. The stellar M/L ratios arc esti- 
mated using a Salpeter (1955) IMF, and using ob- 
served galaxy colors to place constraints on the 
range of plausible stellar M/L ratios, allowing for 
generous variations in galaxy stellar population 
age and metallicity. The B-band luminosities and 
M/L ratios were estimated using stellar popula- 
tion model transformations using the PEGASE 
code (see Fioc & Rocca-Volmerange 1997 for an 
earlier version of the code that we use) . The lines 
show robust fits to the trends, where points more 
than 2.5cr from the trends are excluded from the 
fits, (b) Scrsic index n of elliptical galaxies as a 
function of their absolute magnitude. 



6 




Fig. 2. — (a) Estimated virial coefficient C2 for 
a Sersic model based on an aperture radius re/8; 
C2 is obtainted by using the velocity radial profile 
(Tr{r), the projected radial profile (Tp(i?) and the 
projected aperture radial profile (Jap{Rap)- (b) Es- 
timated values of C2 for two different fixed angular 
apertures (l."6 and 2. "2). Overplotted are the ex- 
pected C2-values for different apertures. 
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Fig. 3. — Relation between the M/L ratio and L 
for our selected galaxies. The filled circles repre- 
sent the relation between the M/L ratio and L 
when C2 accounts for variations in the observed 
value of n (i.e. the nonhomology) of elliptical 
galaxies. The open squares represent the relation 
between M/L and L when C2 is to assumed to have 
a constant value (i.e., homology) equal to the value 
used for a King model with rt/rc ~100 proposed 
by BBF. Note that to transform from the units in 
that paper to ours, we must make the transforma- 
tion C2=4.3C2,SBF- 
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